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We present a phase space description of the process of quantum teleportation for a system 
with an A dimensional space of states. For this purpose we define a discrete Wigner function 
which is a minor variation of previously existing ones. This function is useful to represent 
composite quantum system in phase space and to analyze situations where entanglement 
between subsystems is relevant (dimensionality of the space of states of each subsystem is 
arbitrary). We also describe how a direct tomographic measurement of this Wigner function 
can be performed. 



PACS numbers: 03.67.Hk, 03.65.Bz, 03.67.Lx 



I. INTRODUCTION 

Quantum teleportation |Q is a scheme by which an un- 
known quantum state is transported between two parties 
by only transmitting classical information between them. 
This remarkable task requires also that the two parties 
share the two halfs of a composite system prepared in a 
specific entangled state. In recent years this scheme has 
been studied in great detail. Interest in teleportation 
(whose experimental feasibility has also been recently 
demonstrated in simple systems j2|) is motivated by more 
than one reason. In fact, teleportation is one of the most 
remarkable processes requiring the use of entanglement 
as essential resource |]||. Moreover, teleportation can 
also be conceived as a primitive for universal quantum 
computation Q. In the seminal work of Bennett et al 
H the teleportation of the state of a qbit (a two level 
system) was described. Shortly afterwards, Vaidman |^] 
showed how to generalize this method to include systems 
with a space of states of arbitrary dimensionality. In this 
way, the idea of teleporting states of systems with con- 
tinuous variables was first proposed. A concrete analysis 
of the teleportation of the state of a continuous system 
(the electromagnetic field) was first discussed by Braun- 
stein and Kimble ||. In their work, these authors pro- 
posed (and later performed (?]]) interesting experiments 
to accomplish teleportation of continuous variables. For 
this case, it was natural to describe the whole proce- 
dure in terms of phase space distributions ||. However, 
this is not the case for systems with a finite dimensional 
Hilbert space, where the use of the phase space represen- 
tation is not so common. The description of the usual 
teleportation protocol in phase space has been recently 
presented by Koniorczyk, Buzek and Jansky |9) using the 
discrete version of Wigner functions originally introduced 
by Wooters . This approach, as mentioned in [^) , can 
only be used when the dimension of the Hilbert space 
of the system to be teleported is a prime number. In 
this paper we extend the results presented in |9| to the 



case where the space of states has arbitrary dimensional- 
ity. For this purpose we use a different definition for the 
discrete Wigner function that turns out to be very conve- 
nient to analyze situations where entanglement between 
subsystems of arbitrary dimensionality is an important 
issue. 

Several methods exist to represent the quantum state 
of a system with an A dimensional Hilbert space in phase 
space. As mentioned above, Wooters introduced a dis- 
crete version of the Wigner function that has all the de- 
sired properties only when A is a prime number ]l0| . His 
phase space is an A x A grid (if A is prime) and a Carte- 
sian product of such spaces corresponding to prime fac- 
tors of A in the most general case. On the other hand, 
a different approach to define a Wigner function for a 
system with an A dimensional Hilbert space was intro- 
duced by Leonhardt O] that rediscovered results previ- 
ously used by Hannay and Berry Jl3] and others. This 
method, that has the property of being well defined for 



arbitrary values of A, was used in several contexts 14 
and recently applied to analyze the phase space repre- 
sentation of quantum computers and algorithms fllq , |l6[ . 
In this case, the phase space is constructed as a grid of 
2A x 2A points where the state is represented in a redun- 
dant manner (only Ax A of them are truly independent). 
In this paper, we use a hybrid approach allowing us to 
capture the most useful features of both Wooters and 
Leonhardt methods. Thus, to represent a quantum state 
of a bipartite system we use, following Wooters, a phase 
space which is a Cartesian product of two grids. Each 
one of these grids has, as it does in Leonhardt approach, 
2Ni x 2Ni points (where Aj is the dimensionality of the 
Hilbert space of the i-th subsystem) . 

The paper is organized as follows: In Section II we first 
review the usual approach to define Wigner functions for 
systems with an A dimensional Hilbert space. Then 
we introduce a convenient generalization that can be 
adopted in order to analyze bipartite (or multi-partite) 
systems. We discuss some general properties of the 
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Wigner function and analyze the phase space represen- 
tation of a family of entangled states (generalized Bell- 
EPR states). In Section 3 we show how to describe tele- 
portation of the quantum state of a system with an N 
dimensional Hilbert space using the Wigner function. In 
Section IV we describe the procedure by which a direct 
measurement of the Wigner function of a composite sys- 
tem can be done. In Section V we summarize our con- 
clusions. 



II. PHASE SPACE REPRESENTATION OF 
COMPOSITE QUANTUM SYSTEMS 

We describe here the formalism of Wigner functions to 
represent a composite quantum system in phase space. 
For simplicity, we consider the composite system to be 
formed by subsystems each one of which has an N di- 
mensional Hilbert space We first describe the properties 
of discrete Wigner functions for one of the subsystems 
and later discuss the phase space representation of the 
composite system. It is clear that one can always split 
a system into subsystems in many ways. For example, 
when TV is a composite number one can choose to sep- 
arate each subsystem into even smaller subsystems each 
one of which has a space of states with the dimensional- 
ity of the prime factors of N. Adopting this description 
is simply a matter of physical convenience. Here, we as- 
sume that the N dimensional subsystems are the relevant 
elementary components and that, in some physically in- 
teresting regime, the entanglement between them can be 
manipulated. 



A. Discrete Wigner functions 

To represent the quantum state in phase space we must 
first define the notions of position and momentum. To 
do this, we introduce a basis of the Hilbert space B x = 
{|n),n = 0, ...,N— 1} which we arbitrarily interpret 
as the position basis (with periodic boundary conditions: 
\n+N) = \n)). Given the position basis B x , we introduce 
the conjugate momentum basis B p = {\k), k = 0, . . . , N— 
1} by means of the discrete Fourier transform. The states 
of Bp can be obtained from those of B x as 



\k) = -LV exp(i2 



imk/N)\n). 



(1) 



As in the continuous case, position and momentum are 
related by the discrete Fourier transform. The correct 
semi-classical limit corresponds to the large N limit since 
the dimensionality of the Hilbert space is related to an 
effective Planck constant as N = l/2%h. 

Displacement operators in position and momentum, 
denoted as U and V, are defined as: |Tfl ]: 



U m \n) = \n + m), (j m \k) = exv{-2mmk / N)\k) , 
V m \k) = \k + m), V m \n) = exp(i2irmn/N)\n). (2) 

Commutation relations between U and V directly gener- 
alize the ones corresponding to finite translations in the 
continuous case: 



yvfjq = UiVPexp(i2Trpq/N). 



(3) 



A reflection operator R can also be defined as the one 
acting in the position basis as R\n) = \ — n) (again, 
this operation is to be understood mod N). R is re- 
lated to the Fourier transform Uft (where (n'\UFT\n) — 
exp(i2irnn' / N)) since R — Uprp. 

To represent the state in phase space we use the 
Wigner function defined as the following expectation 
value (110: 



W(a) = Tr(i(a)/3) 



(4) 



where a denotes a phase space point (a = (q,p)) and 
A(a) are the so-called "phase space point operators" de- 
fined in terms of displacements and reflections as 116 11 



A(a) = ^U q RV~ p exp(inpq/N). 



(5) 



It is important to mention that, in order for W(a) to have 
all the desired properties the phase space should be a grid 
of 2N x 2N points which for the rest of the paper will be 
denoted as G2N (i- e -j GW is the set of points a = (q,p) 
where q and p take values between and 2N — 1). It will 
also be useful to denote the first N x N sub-grid as Gat 
(i.e., Gat is the set of points a = (q,p) where q and p 
take values between and N — 1). 

The Wigner function (^) obeys three defining proper- 
ties: (PI) It is real valued, which is a consequence of the 
fact that the operators A(a) are hermitian by construc- 
tion. (P2) The Wigner function can be used to compute 
expectation values between states as 



Tr[p APB ]=N W A {a)W B (a) 



(6) 



This is a consequence of the completeness of the set of 
operators A (a), that satisfy 

Tr[A(a)A(a')} = -L 5 N (q> - q)5 N (p> p) (7) 

where a, a' e G N and S N {q) = J2n=o e 2 ™ 9 ^ is 
the periodic delta function (which is zero unless q = 
mod N). 

As A(a) form a complete orthonormal basis of the 
space of operators, one can expand the density matrix 
p in this basis and show that the Wigner function are 
nothing but the coefficients of such expansion. Thus: 
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p = 4N W{a)A(a) (8) 
= N W(a)A(a). (9) 

The last expression, where the sum is over a 6 GW, can 
be obtained from (|8j) by noticing that the contribution of 
each of the four N x N sub-grids are identical. This can 
be shown by using the fact that (for a q , a p — 0, 1): 

A(q + <7 q N lP + a p N) = A(q,p) (-lypl+W+WpN . 

(10) 

Finally, H^a) satisfies a third crucial property. Con- 
sider a line L in the grid G 2 n (a line L is the set of all 
points a — (g, p) such that ap — bq = c for given inte- 
gers a, b and c). Then the sum of W^(a) over all points 
lying on any line L is always positive. This property 
(P3) is a consequence of the following fact: Adding all 
phase space operators over a line L defined by the equa- 
tion ap — bq = c one obtains a projection operator. Thus 
Al = Y^aeL A(a) is a projector onto an eigenspace of 
the phase space translation operator 

f(a, b) = U a V b exp(i7ra6/iV) (11) 

with eigenvalue exp(inc/N). (If T(a,b) does not have 
exp(iirc/N) as one of its eigenvalues the projector Al 
vanishes). The simplest example of the use of this prop- 
erty is the following: consider the horizontal lines p = c. 
The sum of W^a) over these lines vanishes if c is odd 
(because T(1,0) = U has eigenvalues exp(ink/N) when 
k is an even integer). On the other hand, when c is 
even, the sum of W(a) is equal to the probability of mea- 
suring a momentum equal to c/2. Thus, ^2 q W(q,p) = 
(p/2\p\p/2) Hp is even (and zero otherwise). 

Let us summarize the results presented so far: The 
Wigner function for systems with an N dimensional 
Hilbert space is defined in (0) as the expectation value of 
the phase space operator A{a) given in ([5]). This defini- 
tion is such that VF(a) is real, it can be used to compute 
inner products between states and it gives all the correct 
marginal distributions when added over any line in the 
phase space, which is a grid G2N with 2N x 2N points. 
The size of the phase space grid is important to obtain 
a Wigner function with all the desired properties. The 
values of W(a) on the sub-grid Gn are enough to re- 
construct the rest of the phase space (since the set A(a) 
is complete when a belong to the grid Gn). However, 
the redundancy introduced by the doubling of the num- 
ber of sites in q and p is essential when one imposes the 
condition that all the marginal distributions should be 
obtained from the Wigner function. 

As an example, we show here the Wigner function of a 
position eigenstate p = |<?o)(9o| (see |16| for other inter- 
esting examples): 

W(q,p) = ^-6 N (q - 2q ) (-l)P<«-=*>". (12) 



Here, Zn denotes z modulo N. W(a) is non-zero only on 
two vertical strips located at q — 2q modulo N. When 
q = 2q we have W(2q Q ,p) = 1/2N while for q = 2q ±N 
the modulating factor (— l) p produces oscillations. They 
can be interpreted as the interference between the q = 2qo 
strip and its mirror image induced by the periodic bound- 
ary conditions. The fact that W(a) becomes negative 
in this interference strip is essential to recover the cor- 
rect marginal distributions. Adding the values of W(q,p) 
along a vertical line gives the probability of measuring 
q/2 which should be 1 for q — 2q Q and zero otherwise. 

B. Wigner function for composite systems 

We now consider a composite system with Hilbert 
space HS 1,2 ^ — TiS 1 ' ® HS 1 ^ (for simplicity we assume 
that the dimension of both spaces US' 1 ' is the same but 
the method can be generalized if this is not the case). 
To represent the states of this composite system in phase 
space we generalize the approach described in the previ- 
ous subsection in an obvious way: We consider the phase 
space for the composite system as the Cartesian product 
of the ones for the subsystems (as in the classical case) 
and use the phase space point operators defined as 

A{a u a 2 ) = A(ai) ®A(a 2 ). (13) 

The set {A(ai, a 2 ), with 014 £ Gn}, is a complete or- 
thonormal basis of the space of operators on Ti 1 - 1 ' 2 ^ since: 

Tr h2 (A(a 1 ,a 2 )A(a' 1 ,a 2 )) = — — ^ fofai - a[) 

x 5 N (a 2 - a' 2 ), (14) 

where ai,a^ € Gn and Sn{o) = Sn(q)Sn(p)- 

The Wigner function for the composite system is de- 
fined as the expectation value of the above operators: 

W(a 1 ,a 2 ) = Tr(A(a u a 2 )p) (15) 

This function has the same properties than the one for 
each subsystem. In fact, the three properties (P1-P3) 
generalize trivially to this case. Reality (PI) is once 
again an obvious consequence of the hermitian nature 
of the phase space point operators. Completeness of the 
operators A(a\, a 2 ) enable us to expand the total density 
matrix in this basis and write: 

P={AN) 2 W(a u a 2 )A(a 1 ,a 2 ), 

cxi ,0:2 £Gn 

= N 2 W(ai,a 2 )A(ai,a 2 ). (16) 

Oil ,CK2 f=G27V 

Once again, the first line (where both a\,a 2 € Gn) can 
be transformed into the second line where the sum can 
be extended to the grid G 2 n by noticing that the contri- 
bution over sub-grids are identical (due to the relations 
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obtained by using (|i~0|)). Using this, it is simple to show 
that inner products between two states of the composite 
system can be computed from the Wigner functions as 
(P2): 

Tr(pAPB) = N 2 W A {a x ,a 2 )W B {a x ,a 2 ). (17) 

The third property (P3) is valid as well but it is worth 
discussing it with more detail. The Wigner function 
turns out to be positive when added over any "slice" of 
the total phase space. A "slice" in phase space can be 
defined (following Wooters jlO]]) as a set of lines {Li, L 2 } 
(one line for each subsysteniJTThese sets are called slices 
since in the continuous limit, the set of phase space points 
satisfying the equations aip\ — biqi = c\ (that defines L{) 
and a 2 p 2 — b 2 q 2 = C2 (defining L 2 ) form a two dimensional 
manifold immersed in the four dimensional phase space. 
The fact that the Wigner function is positive when added 
over all points belonging to the slice {Li,L 2 } is obvious: 
Thus, adding A(a\,a 2 ) over all points where ct\ and a 2 
respectively belong to the lines L\ and L 2 one always ob- 
tains a tensor product of two projectors. These operators 
project onto the eigenspaces of the operators 

T 1 =T(a 1 ,b 1 )®I, 
T 2 =l®T{a 2 ,b 2 ) 

with eigenvalue exp(i7rci/Af) and exjp(iirc 2 /N) respec- 
tively. We omit the proof of this property, that can be 
done by copying the one presented in jl6| (see also below 
for a simple proof of a related property). Slices defined 
as above, by picking a line for each subsystem, will be de- 
noted "separable slices" . As we have just seen, Wigner 
functions when added over separable slices can be used to 
compute probabilities for the outcomes of measurements 
of separable observables. These are properties that are 
measured by independent experiments on the two sepa- 
rate subsystems. 

It is possible, however, to consider more general non- 
separable slices in the phase space. As we will show 
below, by adding phase space operators on these non- 
separable slices we will obtain projectors onto entangled 
states. For this we should use a general kind of mani- 
fold Li j2 that can be defined as the set of points (ai,a 2 ) 
satisfying the condition a±pi — b\q\ + a 2 p 2 — b 2 q 2 = c\ 2 
(notice that calling L\ 2 a line can be misleading: in the 
continuous limit it is a three dimensional manifold). It 
is simple to show that by adding phase space operators 
over all points belonging to L\ t2 one also obtains a projec- 
tion operator. It is rather instructive to see this proof: 
For this, we just have to use the fact that the Fourier 
transform of the subsystem's phase space point opera- 
tors A(oti) is a translation Jig : 

2JV-1 2jT 
T(a,b) = MQ>P) exp(-i— (ap-6g)). (18) 

q,p=0 



Using this, we can compute the sum of A(ai,a 2 ) over 
Li 2 as follows: 

(01,02)6^1,2 

Y i(ai) <g> A(a 2 ) 

xS 2N (aipi - hqi + a 2 p 2 - b 2 q 2 - c 12 ) 

2N-1 

= W E T A (a 1 ,6 1 )®r A (a 2 ,6 2 )e^ Aci2 , (19) 

where to obtain the last line we represented the delta- 
function as a sum of exponentials and used (pf). As 
the translation operators are unitary and cyclic one can 
always express them in terms of their eigenstates and 
eigenvalues as 

N-l 

T(a t ,b t ) = \</>i,n)(<t>i,n\eM-tt™/N). (20) 

n=0 

Using this equation one finally obtains: 

Ali,2 = \<f>l,n) {<f>l,n\ ® \(j>2,m){<j>2,m\ 

X 5 2N (n + m- ci )2 /2). (21) 

This explicitly shows that Al 1 2 is a projector onto the 
eigenspaces of a collective operator of the bipartite sys- 
tem. Thus, ( pi] ) projects onto eigenspaces of 

T , i, 2 = f(ai,6i)®f(a 2 ,62) (22) 

with eigenvalue exp(i7rci,2/A r ). It is clear that, when N 
is even the projector is non-vanishing only if ci,2 is even. 
The most important conclusion is that, generically, the 
sum of phase space point operators over non-separable 
manifold Li 2 will correspond to a projector over an en- 
tangled state. In the previous case T\ >2 is a collective op- 
erator generating simultaneous phase space translations 
of both systems. In the continuous limit it is clear that 
this operator is generated by linear combinations of the 
momenta and coordinates of the two subsystems. 

In general, a manifold L\ 2 will not be associated with 
a one dimensional projector. For example, the manifold 
L +jP/3 defined by the equation p\+p 2 = 2pp corresponds 
to an N dimensional subspace. The same is true for 
the manifold L— )W defined as the set of points satisfy- 
ing qi — q 2 = 2qp. The intersection between these two 
sets will be denoted as Lp (defined as the set of points 
belonging to both £+, P/3 and L- } q„) and corresponds to 
a one dimensional projector over an entangled state (see 
below) . 

Wigner functions of separable states have very dif- 
ferent features than those of entangled states: In fact, 
if the density matrix is a tensor product p — p\ ® p 2 
then W{a\,a 2 ) is a product of the form W{ai,a 2 ) = 
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W\{a\)W2{ci2) ■ More generally, if the state is separa- 
ble (i.e., the density matrix is a convex sum of tensor 
products) then the Wigner function is a convex sum of 
products of the above form. For entangled states (states 
which are not separable) this is not the case as will be ex- 
plicitly seen below. It is also useful to notice that reduced 
Wigner functions can be computed for one subsystem by 
adding the total Wigner function over the complemen- 
tary half of the phase space. Thus, summing a 2 ) 
over a2 is equivalent to taking the partial trace over the 
second subsystem since 

Wi(ai)= W(ai.aa) 

02EG2N 

= Tr 1 (i(a 1 )p 1 ), (23) 

where p\ is the reduced density matrix of the first system 
obtained from the total density matrix as p\ = Tr 2 p. 
Finally, it is also useful to notice other properties of the 
composite Wigner function when we sum it over half of 
the phase space. For example, if pa and ps are two states 
of the composite system, then 

F(a 1 ,a' 1 )= ^2 W A {ai,a 2 )WB(a' 1 ,a2) 

= i ^Tr 2 (Tr 1 (p A i(a 1 ))Tr 1 (p s i(a' 1 )), (24) 
an equation that will useful later. 

C. Wigner functions for Bell states 

Let us first introduce a complete basis of entangled 
states (generalized Bell states Q). First we define the 
state | Go) as 

i 

l °> = -/^El n >i®l n >2- ( 25 ) 

v n— 

This pure state for the composite system is maximally 
entangled since the reduced density matrix of each sub- 
system is proportional to the identity. A complete basis 
of entangled states can be defined from |9o) as follows: 

1 N ~ 1 

= -7= E ^ P " n,N \n)i ® \n - qp)2. (26) 

v n— 

These states can all be obtained from |Oq) by using one 
of the following equivalent expression: 

\Qp) =T l {q^pp)^^' N ®l 2 \Q Q ) 

= v?*®ar w |eo) 

(27) 

Above, the notation (3 — (q@,pp) was used. When j3 £ 
Gn these states form a complete orthonormal basis of the 
Hilbcrt space HS 1,2 ^ . In fact, these N 2 vectors satisfy 



(QplQp) =5 N (/3-l3'). (28) 

It is interesting to notice that these are the common 
eigenstates of the following complete set of commuting 
operators: 

U+ = U t ® U 2 , 

V- = V X ®V%. (29) 

In fact, V- and U+ commute and that Bell states ( |26| ) 
obey the following identities: 

U+\Q P ) = exp(-i2np fi /N)\Q p ), 

V-\O ) = exp(i27rq /N)\Q ). (30) 

These expressions allow us to better understand the 
nature of Bell states: U+ displaces both systems in po- 
sition by the same amount while V- displaces them in 
momentum in opposite direction. As Bell states are 
common eigenstates of these operators, such states can 
be interpreted as corresponding to the eigenstates of 
the toted momentum and relative position (note that in 
the continuum limit [/+ = cxp(— i(P± + P2)8x/K) and 
V- = exp(iSp(Qi — Q2)/K). In this sense, Bell states 
( p6| ) are precisely the ones used by Einstein, Podolsky 
and Rosen |isj | to present their argument against com- 
pleteness of quantum mechanics as a description of na- 
ture. 

Having this in mind, one expects the phase space rep- 
resentation of Bell states to be rather simple. This is 
indeed the case: 

W p {a u a 2 ) = Tr(\Q ){Q \ A( ai ) <g> A(a 2 )) 
= Wo(a 1 -2j9, a 2 ) 

where 

W {a) =Tr(|6 )(e |i(a 1 )®i(a 2 )) 

= (2jV) 2 < ^ gai ~ qa ^ SN ( Pai + P a ^ 

Thus, the Wigner function of \Qp) is nonzero only when 
qa 1 —qa 2 = 2(?/3 and p ai +p a2 — %Pf3 (modulo AT). The os- 
cillations, whose origin we described above for a simpler 
case (|l2l), can also be interpreted as due to the inter- 
ference with the mirror images created by the boundary 
conditions . Notice that these are precisely the equations 
defining the manifold Lp. Thus, the projectors onto Bell 
state 1 6/3) is simply the sum of phase space point opera- 
tors over all points belonging to Lp. 

Before explicitly discussing teleportation it is useful 
to present some further results related to Bell states 
and their connection to phase space point operators. A 
complete basis of the space of operators on H^ 12 ^ can 
be constructed using Bell states: Thus, the operators 
S(/3i, /3 2 ) = |0/3i)(©/3 2 | form a complete orthogonal set 
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(with /3i, fa S Gn)- The change of basis between this set 
and phase space point operators is 

B(0i,0 2 ) = (ANf K{0i,0 2 \ai,a 2 )A(ai,a 2 ). 

(32) 

where the coefficients K(0i, 2 \cti, ct 2 ) are, in gen- 
eral, complex numbers (when 0\ = 2 we have 
K{0,0\ai,a 2 ) — Wp(cti,a 2 )). The precise form of this 
coefficients can be easily obtained but will not be needed 
here. The inverse of (|3^) is also useful and reads 

A(ai,a 2 )= K{ai,a 2 \0i,0 2 )B(0i,0 2 ) (33) 
Pu Pi 

Simple relations between the coefficients of ([32]) and d33| ) 
exist. In particular, one can show that 

K(ai,a 2 \0i,0 2 ) = K(0 2 ,0i\ai,a 2 ) 

= K*[0i,0 2 \ai,a 2 ) (34) 

Finally, we mention that Bell states satisfy the follow- 
ing identity: 

Tr 1 (i(a)Ol 2 |e )(e |) = ^l 1 ®A T {a), (35) 

where the transpose of the phase space point operator 
appears in the right hand side. 

III. TELEPORTATION IN PHASE SPACE 

We show here how the usual teleportation protocol can 
be described in phase space. We consider three identical 
subsystems labeled by the integers j = 1,2,3. The aim 
is to teleport the initial state of system 1, which is char- 
acterized by an arbitrary Wigner function W(ai). For 
this we initially prepare systems 2 and 3 in one of the 
Bell states (for simplicity we use 1 60)2.3 as the initial 
state). Thus, the initial density matrix of the combined 
three-partite system is 

Pi, 2, 3 = Pi ® |@o)(@o|2,3 

= N 3 W(ai)W {a 2 ,a 3 ) 

Q'l ,Q2 -OtZ^G2N 

xl(ai) ®A{a 2 ) <g)i(a 3 ) (36) 

where the Wigner function for the Bell state 10^) is given 
in (|3l|). After preparing this initial state the teleporta- 
tion protocol proceeds as follows: First we perform a 
measurement of system 1 and 2 that projects them into 
the Bell basis. Physically, as discussed above, this corre- 
sponds to a collective measurement that determines the 
total momentum pi+p 2 and the relative coordinate qi—q 2 
for these two subsystems. After this measurement the 
state of the system is projected into the resulting state 



|0/3)i,2, where = (qpiPp) {pp and qp are the measured 
values of the total momentum and distance). To explic- 
itly write down the resulting state, it is better to rewrite 
equation ( j36"l) expressing the phase space point operator 
A(cti) (g> A(a 2 ) in terms of the Bell operators B(0i,0 2 ) 
as in (||). Thus, 

Pi.2,3 = (4A0 3 W(ai)W o (a 2 ,a 3 )K(ai,a 2 \0i,0 2 ) 

xB(0i,0 2 )®A(a 3 ) (37) 

From this equation it is obvious that, after the Bell mea- 
surement of the first two subsystems, only the terms with 
01 = 2 = in the above expression survive. Therefore, 
the state of the third subsystem is 

p' 3 =4N ^'(03)^(03) (38) 

a 3 SGjv 

where the term multiplying the phase space point opera- 
tor in ( [37| ) can be identified as the new Wigner function 
of the third system (up to a normalization constant). In 
turn, this Wigner function can be seen to be defined by 
an expression involving a sum that contains the initial 
Wigner function of the first system. Thus, this can be 
written as 

W'(a 3 )= J2 Z{a z , ai )W{ ai ). (39) 

aiEGjv 

The matrix Z(a 3 , a\) simply tells us how to construct the 
final Wigner function for the third system from the initial 
Wigner function of the first one. Again, the explicit form 
of this matrix is read from ( |37T) : 

Z(a 3 ,ai) = (2N) 4 ^ W (a a , a 2 )Wp{ai, a 2 ). (40) 

The above expression is easy to evaluate since it contains 
a sum over half of the total phase space. Thus, we can 
use the relation ( p4] ) to simplify it. Doing this (and using 
(p5|)) one discovers that Z(a3,ai) is just the trace of a 
product of two phase space point operators acting on the 
second subsystem which are evaluated at points ot\ and 
as~ 20- Therefore, taking into account the orthogonality 
conditions for phase space point operators, one obtains 

Z(a 3l ai) = S N {a 3 - ai - 20). (41) 

This means that the state of the third subsystem has a 
Wigner function that is displaced in phase space by an 
amount with respect to the initial state of the first 
system. Therefore, to recover the initial state one must 
displace the third system in phase space by applying the 
evolution operator T{0) — T(qp,p/3) (in fact, one can 
show that when the operator T{0) acts on the state of 
the system, the Wigner function is simply shifted in phase 
space by the amounts (2a, 26) the factors of 2 being orig- 
inated in the fact that the grid has 2N x 2N points jl6| ) . 
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Obviously, the recovery operation depends on [3, which is 
revealed by the Bell measurement performed on the first 
two subsystems. In this way the final state of the third 
system is always identical to the unknown initial state of 
the first subsystem, which is the goal of the teleportation 
protocol. 



IV. MEASURING THE WIGNER FUNCTION 

The Wigner function, as mentioned above, is in one 
to one correspondence with the quantum state of the 
system. Thus, experimentally determining the value of 
W(ai,ot2) in every phase space point is equivalent to 
completely determining the state of the system. Exper- 
imentally reconstructing the Wigner function has been 
the goal of a series of experiments, all dealing with con- 
tinuous systems JTgj ] . In general, these experiments are 
aimed at determining first marginal distribution of some 
observables and later reconstructing the Wigner function 
by means of a Radon-like transform. However, it is in- 
teresting (and useful) to realize that there is no need to 
fully determine the quantum state to evaluate the Wigner 
function at a given phase space point. Indeed, this was 
proposed originally by Davidovich and Luttcrbach for the 
Wigner function of the state of the electromagnetic field 
stored in a superconducting cavity [KJ. Using a minor 
variation of this method the direct measurement of the 
Wigner function at the origin of phase space has been 
recently performed in a cavity QED experiment. More 
recently, this method was generalized and shown to be 
applicable to the measurement of discrete (or continuous) 
Wigner functions of generic systems in Jl6| , p2[ . Here, we 
will describe how this tomographic scheme can be gen- 
eralized to directly measure the Wigner function we pre- 
sented in this paper. The fact that this Wigner function 
can be measured in an efficient way (i.e., in a number 
of steps which scales polynomially with the dimension- 
ality of the system) may be important if, for example, 
one is able to relate interesting physical quantities (like 
some entanglement measures) with phse space observ- 
ables (work is in progress in this direction). 

The efficient strategy to measure the Wigner func- 
tion of a composite system at (any) given phase space 
point is, as mentioned above, a direct generalization of 
the idea originally proposed in jl6],^2| to measure the 
Wigner function of an N dimensional system. The ba- 
sic ingredient can be described in terms of the following 
quantum algorithm. Consider a system initially prepared 
in a quantum state p. We put this system in contact 
with an ancillary qbit prepared in the state |0). This 
ancillary qbit plays the role of a "probe particle" in a 
scattering-like experiment. The algorithm is: i) Ap- 
ply an Hadamard transform to the ancillary qbit (where 
2T|0> = QO) + |1))/V2, H\l) = (|0) - |1))/V2), ii) Apply 
a "controlled-M" operator (if the ancilla is in state |0) 
this operator acts as the identity for the system but if the 



state of the ancilla is |1) it acts as the unitary operator 
M on the system), iii) Apply another Hadamard gate to 
the ancilla and finally perform a weak measurement on 
this qbit detecting its polarization (i.e., measuring the 
expectation values of Pauli operators a z and <j y ). It is 
easy to show that the above algorithm has the following 
remarkable property: 

(<r z ) = Re(Tr(Mp)), (a y ) = Im(Tr(Mp)). (42) 

Thus, the final polarization measurement of the ancillary 
qbit reveals a property determined both by the initial 
state p and the unitary operator M . 

In |22| we discussed how to view this simple algorithm 
as the basic tool to construct a rather general tomogra- 
pher (and also a rather general spectrometer). In par- 
ticular, we showed how to use it to measure the Wigner 
function of a simple system. Here, we show how to adapt 
it to measure the Wigner function of the composite sys- 
tem we have been discussing so far. This can be done by 
applying the algorithm shown in Figure 1. 



10x01 
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Measure <o, > 



FIG. 1. Circuit for measuring W (qi, 02), for a compos- 
ite system. The controlled operations apply the operator 
2NA(ai) (which is unitary) conditioned on the state of the an- 
cilla qbit. The phase space operators parametrically depend 
on the phase space point on and can be efficiently built as a 
simple sequence of displacements and reflections on each of 
the two subsystems. The measured polarization of the ancilla 
qbit is directly proportional to the Wigner function. 

From the previous discussion it is clear that the above 
algorithm is such that by measuring the polarization of 
the ancillary qbit we determine the Wigner function. In- 
deed, this follows from the identity 

(a z ) = AN 2 Tv{pA{ai) <g> A{a 2 )) = AN 2 W(a l ,a 2 ). (43) 

As phase space point operators ([5]) are simply a prod- 
uct of displacement operators (which implement addition 
of one, modulo N) and reflections (which are the square 
of the Fourier transform) the network of Figure I can 
be implemented efficiently (i.e., it involves a number of 
elementary gates that grows polynomially with log (A)). 
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V. CONCLUSION 

In this paper we used a hybrid approach to construct 
a Wigner function to represent quantum states of a com- 
posite system in phase space. The function we defined 
has interesting features enabling us to study situations 
where entanglement between subsystems plays an impor- 
tant role. Thus, the hybrid method captures some of the 
most useful properties of the Wigner functions defined 
by Wooters Jin] and Leonhardt (and others) For 
a bipartite system this function depends upon two phase 
space coordinates W(ai, 0:2 )■ The phase space is a Carte- 
sian product of the phase spaces of the subsystems, as it is 
the case in Wotters proposal. However, each phase space 
grid has 2N x 2N points, as suggested by Leonhardt and 
others. For separable states the Wigner function is, in 
general, a convex sum of products of independent func- 
tions for each subsystem. Thus, this Wigner function 
is a natural tool to study entanglement between subsys- 
tems. In this paper we showed that basis of entangled 
states can be identified with non-separable slices in the 
phase space (the basis formed by Bell states is one such 
example). We also showed that W(a, a') is measurable 
by a simple scattering-like experiment where an ancillary 
particle successively interacts with the two subsystems. 
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